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We present the results of <j> meson production in the K + K~ decay channel from Au+Au collisions 
at y / ijvjv = 200 GeV as measured at mid-rapidity by the PHENIX detector at RHIC. Precision 
resonance centroid and width values are extracted as a function of collision centrality. No significant 
variation from the PDG accepted values is observed, contrary to some model predictions. The <j> 
transverse mass spectra are fitted with a linear exponential function for which the derived inverse 
slope parameter is seen to be constant as a function of centrality. However, when these data are fitted 
by a hydrodynamic model the result is that the centrality-dependent freeze-out temperature and 
the expansion velocity values are consistent with the values previously derived from fitting identified 
charged hadron data. As a function of transverse momentum the collisions scaled peripheral-to- 
central yield ratio Rcp for the <j> is comparable to that of pions rather than that of protons. This 
result lends support to theoretical models which distinguish between baryons and mesons instead 
of particle mass for explaining the anomalous (anti)proton yield. 
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I. INTRODUCTION 

Relativistic heavy-ion experiments have a goal of pro- 
ducing matter at extreme temperatures and energy densi- 
ties such that conditions are favorable for the transition 
to a deconfined state of quarks and gluons, the Quark 
Gluon Plasma (QGP). Theoretical calculations predict 
that the temperatures and energy densities which can be 
reached at the Brookhaven National Laboratory's Rel- 
ativistic Heavy Ion Collider (RHIC) will exceed those 
needed for the formation of the QGP 0, HE1 • 

The production and decay of the (f> meson have long 
been recognized as an important probe for the state of 
matter produced in relativistic heavy ion (RHI) colli- 

^rUJJI S S S S M M III El III EE Ellll 

1H l2l u3 • In PP collisions the creation of the 4> is sup- 
pressed according to the Okubo-Zweig-Iizuka rule |22fl . 
Hence, if there is an enhancement of the </> yield in 
RHI collisions relative to pp collisions, this could be ev- 
idence of non-conventional production mechanisms such 
as strange quark coalescence via the formation and sub- 
sequent hadronization of the QGP. The fact that the 4> 
yield is undistorted by feed-down from higher mass reso- 
nances makes it an attractive probe in this respect. 

The decay modes of the <fi, specifically the dilepton 
channels (e + e~ or and the K + K~ channel, will 

probe the final state differently should the decay take 
place in the presence of the QGP-mixed or the com- 
pletely hadronized phase. The dileptons will have in- 
significant interactions with the medium, while the kaons 
can scatter until freeze-out. The lifetime of the 4> in vac- 
uum is large (~ 45 fm/c) compared to say a 10 fm sized 
interaction region. However, several theoretical calcula- 
tions [T3L I2TL l23| predict that the c/> mass and width could 
be significantly modified in either the hot or the cold nu- 
clear medium. These medium induced effects could be 
manifested through measured shifts in the mass centroid 
of the resonance or changes in the resonance width. Also 
predicted are changes in the relative branching ratio be- 
tween kaon and lepton pairs, with respect to the Particle 
Data Group (PDG |3). 

The production mechanism of strangeness in heavy ion 
collisions can be investigated through the measurement 
of the particle yields. In this paper, we study system 
size dependence by analyzing centrality selected data. A 
comparison between different systems can be made by 
normalizing to the number of participant pairs. The ex- 
pectation is that for production dominated by soft pro- 
cesses, the yields scale as the number of participants. 
We compare the centrality dependence of strange and 
non-strange particle yields in order to reveal the possible 
flavor dependence. 

An additional important question is whether the (\> 
mesons participate in radial flow together with the other 
hadrons, or if they freeze-out earlier, as might be true if 



the small vacuum cross sections of the <f> with hadrons 
persist in the fireball. Previous measurements have 
yielded contradictory results [H |2(|. One of the im- 
portant advantages of RHIC experiments is the capabil- 
ity to examine the momentum spectrum as a function of 
centrality for a variety of hadrons which should yield im- 
portant additional information on the radial flow issue. 
A spectral shape analysis including a simultaneous treat- 
ment of the 4> and the more abundant hadrons (jr,K,p) 
will be presented here. 

At high pt, hadrons are primarily produced from the 
fragmentation of hard-scattered partons. One of the 
most exciting results from RHIC was the discovery of 
hadron suppression in central Au+Au collisions I27L l2p 
where this suppression is absent in d+Au collisions |29| . 
Surprisingly, it was also discovered that proton and anti- 
proton production at intermediate px (1.5-4.5 GeV/c) 
scales with the number of binary nucleon-nucleon colli- 
sions (N co u) as would be expected for hard-scattering in 
the absence of any nuclear modification [3(j. In fact, the 
intermediate px anti-proton to pion ratios were found to 
exceed by a factor of 3 the values expected from par- 
ton fragmentation |30l l3l| . These experimental results 
lead to the conclusion that protons and pions have dif- 
ferent production mechanisms at intermediate px |30| . 
One possible explanation invokes parton recombination 
from the QGP [H |H H HI ■ A measurement of the 
nuclear modification factor for the (j) meson, which has 
a mass comparable to the proton but carries only two 
quarks, is crucial for understanding the hadron produc- 
tion at intermediate px- In this work we have measured 
the nuclear modification factor through the ratio Rcp 
of central to peripheral yields scaled by their respective 
N co a value. 

To put our results into perspective, we begin by de- 
scribing the currently available <fi data obtained in RHI 
collisions. The production of <fi mesons has been stud- 
ied systematically at ever increasing ^/s from the AGS 
to RHIC. The E802 collaboration made the first ob- 
servation of the 4> m fixed target central collisions of 
U.6A GeV Si+Au (y^Tjv = 5.39 GeV) via the K+R- 
channel [3(|. They obtained a ratio N^/Nr- = 11.6%, 
roughly consistent with the ratio obtained in pp data over 
a wide range of a/s [j^i . The analysis of the E802 rapidity 
distributions indicated that the </> production scaled with 
the product of the K + and K~ separate yields, and that 
there was either significant rescattering of the (f> after pro- 
duction or the production itself came after rescattering 
of the colliding participants. 

Also at the AGS the E917 experiment has re- 
ported another (f> measurement with 11. 7A GeV Au+Au 
(y/SNN = 4.87 GeV) in the rapidity range 1.2 < y < 1.6 
in five centrality bins 38] The observed yield of the <fi in- 
creased towards more central collisions with a distinctly 
faster than linear dependence on the number of partic- 
ipants. The yield increase of the <fi in central collisions 
was stronger than that of the ir since the N$ /N„ ratio in- 
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creased in central collisions. However, the N^/N K + and 
N^/N K - ratios were constant as a function of centrality. 

At the SPS the NA49 experiment has measured <f> 
production in pp, p+Pb, and Pb+Pb collisions with 
E beam =158A GeV (^/iww = 17.5 GeV) in the rapid- 
ity range 3.0 < y < 3.8 |33. Relative to the pp yields, 
these data showed that the ratio of the <j> yield to the 7r 
yield in central Pb+Pb collisions was enhanced by a fac- 
tor 3.0 ± 0.7. Another SPS collaboration, NA38/NA50, 
has measured the cf> in the fi + fi~ channel El), for 
which the extracted effective temperature and dN / dy dif- 
fer from those obtained in the K + K ~ channel in the same 
systems. The yield difference between the two SPS ex- 
periments has been calculated to be factors of 2-4 |4^ . 
with the NA38/NA50 result being higher. 

The first measurement of the <fi meson at RHIC was 
reported by the STAR collaboration [43 , in the collisions 
of Au+Au at y / s]vjv=130 GeV at three centralities, 0- 
11%, 11-26%, and 26-85% in the rapidity range -0.5 < 
y < +0.5. The extracted temperature T and the ratio 
Nj, /Nh did not vary with centrality. 

One may summarize the current state of knowledge of 
4> production in heavy ion collisions by stating that the 
topic remains highly unexplored territory. The heavy- 
ion measurements do indicate that the observed <f> yield 
is not a simple linear superposition of nucleon-nucleon 
collisions. Rather the data imply the influence of some 
collective effects. Whether those effects are induced by 
cold and/or hot nuclear matter, there do not yet exist 
definitive measurements. Moreover, there has not been 
so far precise enough heavy-ion data which can address 
the question of the change in the <p mass or its width in 
the cold nuclear or the hot QGP medium. And except 
for one experiment measuring the dimuon channel, there 
is a scarcity of useful quantitative information in heavy 
ion collisions concerning the <f> decay into dileptons. 

In this paper we report on a measurement of the <fi yield 
at mid-rapidity in collisions of Au+Au beams from RHIC 
at ^/sjvat = 200 GeV as measured in the K + K~ chan- 
nel by the PHENIX detector. The paper is organized 
as follows. In Section II we give a short review of the 
PHENIX detector configuration. In Section III we de- 
scribe the data analysis procedure. Section IV presents 
and discusses the results. The first precision measure- 
ments of <fr mass and width values obtained in relativistic 
heavy ion collisions as a function of centrality are given 
in Section IV-A. Absolute yields as a function of pr for 
three centrality bins are shown in Section IV-B. The cen- 
trality dependence of yields and ratios are studied in Sec- 
tion IV-D. In Section IV-D the spectra shapes are inter- 
preted in the framework of a hydrodynamical model and 
the freeze-out conditions are extracted. Finally, the nu- 
clear modification factor Rcp for the <fr is obtained and 
compared to those of pions and protons in Section IV-E. 



II. PHENIX DETECTOR 

The PHENIX detector consists of two spectrome- 
ter arms at near zero rapidity, two forward rapidity 
muon spectrometers, and three global event character- 
ization detectors. The central arm spectrometers, shown 
schematically in Fig. are located East and West of 
the beam line with 7r/2 radian azimuthal coverage each. 
These spectrometers are designed to detect photons, elec- 
trons, and charged hadrons. The <f> data for this pa- 
per were obtained with the central arm detector sub- 
systems which provide high resolution particle identifi- 
cation and momentum reconstruction. A complete de- 
scription of the PHENIX appar atus has been published 
elsewhere HI, 0, HI, . We present a brief re- 
view of the relevant detector subsystems in the following 
sections. 



A. Global Detectors 

The global detectors furnish the start time signals, 
collision vertex measurements, and interaction central- 
ity. The centrality for events in the Au+Au collisions 
is determined by combining the data from two sub- 
systems: the Zero Degree Calorimeters (ZDC) and 
the Beam-Beam Counters (BBC) 0. The ZDC are 
hadronic calorimeters located 18 m downstream and up- 
stream of the interaction point along the beam line. 
These calorimeters detect the energy carried by specta- 
tor neutrons. The BBC are Cerenkov telescopes placed 
±1.44 m the center of the beam collision region in the 
pseudo-rapidity region 3.0 < < 3.9. The correla- 
tion between the ZDC energy sum and the charge sum 
recorded in the BBC determines the centrality of the col- 
lision event. The BBC data also determine the longitudi- 
nal collision coordinate (^vertex) an( i the start time for 
the time-of-flight measurements. 



B. Central Arm Detectors 

The central arm spectrometers 0) m PHENIX pro- 
vide charged particle tracking and particle identification. 
This analysis was done with the East arm spectrometer. 
The <j> — > K + K~ decay kinematics are such that that 
the PHENIX detector has negligible acceptance for the 
very low pr 4> particles which would decay into East- West 
kaon pairs. The data included information from the drift 
chamber (DC), the pad chambers (PCI and PC3), the 
high resolution time-of-flight wall (TOF), and the elec- 
tromagnetic calorimeter (EMCal) lead-scintillator detec- 
tors (PbSc), as depicted in Fig. ^ 

Pattern recognition and tracking of the charged par- 
ticles are accomplished using the DC information by a 
combinatorial Hough transform method [HJ. The DC 
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West Beam View 



FIG. 1: Beam's eye view of the PHENIX central arm detector 
subsystems. 



subsystem is located at an average radial distance of 
2.2 m from the beam line. It is a projective tracking 
detector providing high precision measurements in the 
azimuthal XY plane, which are combined with the XY Z 
single spatial point measurement from the PCI located 
at 2.45 m. These data, together with the BBC ^vertex 
information, are sufficient to determine the track's initial 
momentum vector whose magnitude is obtained with a 
resolution Sp/p ~ 0.7% © 1.0% x p (GeV/c). The first 
term in this expression is due to the multiple scattering 
before the DC and the second term is the angular resolu- 
tion of the DC. Based on identified mass measurements 
obtained with the TOF subsystem, the absolute momen- 
tum scale is known to ±0.7%. 

Tracks obtained from the DC/PCI detectors are pro- 
jected to the PC3, TOF, and PbSc detectors where as- 
sociations can be made. The high resolution TOF sub- 
system provides one set of mass measurements while the 
PbSc detectors provide a geometrically independent set 
of mass measurements. The TOF wall is positioned 
5.06 m from the beam line and consists of 960 scintil- 
lator slats oriented along the azimuthal direction. It is 
designed to cover \rj\ < 0.35 and Acf> = ir/4 in azimuthal 
angle. 

The PbSc detector, covering half of the East arm and 
entire West arm, can also be used for hadron timing mea- 
surements. The present analysis uses the PbSc modules 
in the East Arm which are located at 5.1 m in radial 
distance from the beam line and cover a A(f> = 7r/4 az- 
imuthal range. This detector is constructed as separate 
towers of dimension 5.25 x 5.25 x 37 cm 3 , in an alter- 
nating lead-scintillator sandwich type structure ("shish- 
kebob"), approximately 18 radiation lengths in depth. 
As illustrated in Fig. ^ the TOF and the PbSc sectors 
are completely non-overlapping. 



TABLE I: Average number of participants and collisions in 
Au + Au reaction at RHIC for different centralities deter- 
mined from a Glauber model |29|l . The error associated with 
each number is the systematic error. 



Centrality 
(%) 


< Npart > 


< N coU > 


- 10 


325.2 ± 3.3 


955.4 ± 93.6 


10 - 20 


234.6 ± 4.7 


602.6 ± 59.3 


20 - 40 


140.4 ± 4.9 


296.8 ± 31.1 


10 - 40 


171.8 ± 4.8 


398.7 ± 40.5 


40 - 60 


59.9 ± 3.5 


90.6 ± 11.8 


60 - 92 


14.5 ± 2.5 


14.5 ± 4.0 


40 - 92 


32.0 ± 2.9 


45.2 ± 7.3 


Min. bias 


109.1 ±4.1 


257.8 ± 25.4 



III. DATA ANALYSIS 

In this section, we describe the event and track selec- 
tion, particle identification, the details of K + K~ pair re- 
construction, and the corrections for geometrical accep- 
tance, particle decay in flight, multiple scattering, and 
detector occupancy factors, all of which couple into de- 
riving the <j) meson spectra. 



A. Event Selection 

The events selected for this analysis were based on 
the PHENIX minimum-bias trigger provided by the 
beam-beam counters (BBC) and zero-degree calorime- 
ters (ZDC). As noted previously, the centrality of each 
Au + Au collision event was determined by correlating 
the BBC charge sum and the ZDC total energy [5(|. The 
PHENIX minimum-bias data sample included 92.2^'n % 
of the 6.9 barn Au + Au total inelastic cross section |29| . 
This analysis used 20 million minimum-bias events with 
a vertex position within |z V ertexl < 30 cm. 

To study the centrality dependent physics, we divided 
these events into different centrality bins. For the 4> me- 
son line shape analysis, we used five centrality bins: 0- 
10%, 10-20%, 20-40%, 40-60% and 60-92%. The trans- 
verse mass (mr) spectra were reconstructed for three cen- 
trality bins: 0-10%, 10-40% and 40-92%. These bin divi- 
sions were chosen to have approximately equal statistical 
significance for their respective data points. The central- 
ity of collisions is additionally characterized by the aver- 
age number of participants (< N part >) and the average 
number of binary nucleoli nucleoli collisions (< N co u >). 
These two global quantities, shown in Table [I]as a func- 
tion of centrality, are derived from a Glauber model cal- 
culation [2^ . 
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B. Track Selection 

Only tracks with valid information from the DC and 
the PCI were used for the analysis. These tracks were 
then confirmed by matching the projected and associated 
hit information at the TOF wall for the TOF analysis, 
or at the PC3 and the EMCal for the PbSc analysis. 
The differences between the actual azimuthal and lon- 
gitudinal hit coordinates compared to the projected hit 
coordinates were determined. These tracking coordinate 
residuals were converted to standard deviation residu- 
als by a momentum dependent function which computed 
the expected residual coordinate value. On this basis, 
a 3cr track matching cut was used to accept track asso- 
ciations. Lastly, for the TOF wall, an energy loss cut is 
applied on the analog signal height from the scintillator 
slat. This cut has been described in a previous publica- 
tion jH3- 



C. Particle Identification 

As mentioned earlier, the PHENIX central arm spec- 
trometer utilizes the high resolution TOF wall and PbSc 
modules for hadron mass identification. The kaons in 
the TOF wall were identified via reconstructed momen- 
tum combined with a time of flight measurement. The 
timing resolution of this subsystem is a ~ 115 ps. A 
momentum range of 0.3-2.0 GeV/c was selected in order 
to compute the mass distributions 1 . Fig. shows the 
mass-squared distribution of all tracks passing through 
the TOF module for six different momentum bins. The 
kaons were identified by applying a 2a mass-squared cut, 
which is shown by the shaded region in each plot. 

The particle identification with the electromagnetic 
calorimeter modules is sensitive to the fact that the elec- 
tromagnetic and hadronic interactions produce quite dif- 
ferent patt erns of energy sharing between calorimeter 
towers 52] . As a consequence, the hadron timing prop- 
erties of the PbSc depend on the energy deposited on the 
central tower of the cluster, particle momentum, parti- 
cle type, charge, angle of incident of the track, etc. The 
PbSc hadronic timing response was corrected for these 
effects and we obtained an overall timing resolution of 
a ~ 450 ps which is sufficient to enable a clear ir/K 
separation within 0.3 < p (GeV/c) < 1.0 using a 2a 
mass-squared selection criterion. In Fig. [3] the mass- 
squared distributions are plotted for four different mo- 
mentum slices for all tracks passing through PC3 and 
PbSc. The identified kaons are also shown in the figure 



1 The mass-squared for each track is defined as, 
M 2 = p 2 ( i 2 c 2 /L 2 — 1), where p is the momentum, t 
is the time of flight, c is the speed of light and L is the length 
of the path traversed by the track from vertex to the detector. 




M 2 [GevV] 



FIG. 2: Mass-squared distribution of all selected tracks pass- 
ing through the TOF for six selected momentum bins. The 
identified kaons within 2a mass-squared boundaries at the 
different momentum bins are shown by the shaded region on 
each plot. 

by the 2a width shaded histograms superimposed on the 
M 2 distributions for all tracks in different momentum 
bins. 



D. K + K Invariant Mass spectra and <f> Signal 
Extraction 

All identified K + and K~ tracks in a given event were 
combined to form the invariant pair mass distributions. 
Three different pair combinations were used. These are 

a) Both K+ and K~ identified by the TOF detector 
(TOF-TOF combination), 

b) K+ identified with TOF and K~ identified with 
PbSc (TOF-PbSc combination), and 

c) Both K + and K~ identified with PbSc detector 
(PbSc-PbSc combination). 

We did not use the K + from PbSc and K~ from TOF 
in the b) combination. This is due to the fact that 
the PHENIX central arm geometry, in the presence of 
a 0.8 T-m magnetic field, does not have any acceptance 
for such pairs below an invariant mass of 1.06 GeV/c 2 in 
the TOF-PbSc combination. 

A large combinatorial background is inherent to the 
K + K~ pair invariant mass distribution. The combi- 
natorial background was estimated by an event mixing 
method in which all K + tracks from one event were 
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FIG. 3: Mass-squared distribution of all selected tracks pass- 
ing through the PbSc for four momentum bins. The identified 
kaons within 2a mass-squared boundaries at the different mo- 
mentum bins are shown by the shaded region on each plot. 



combined with K ~~ tracks of ten other events within the 
same centrality and vertex class. The mixed event tech- 
nique reproduces the shape of the unlike sign combina- 
torial background. Finally, the size of the combinato- 
rial background is obtained by normalizing the mixed 
event invariant mass spectra to 2 Ak / J ?V ++ ./V__ where N ++ 
and 7V__ represent the measured yields in K + K + and 
K~ K~ mass distributions respectively. This normaliza- 
tion is derived analytically starting from the assumption 
that the number of K ± tracks per event follows a Pois- 
son distribution. A complete derivation of this is given 
in Appendix A. 

The ability of this event mixing method to reproduce 
correctly the shape of the combinatorial background dis- 
tribution was confirmed by constructing, in a similar way, 
the mixed-event like-sign spectrum and comparing it to 
the same-event like-sign pair distribution. The assump- 
tion is that the like-sign pair distributions are purely 
combinatoric. For the three detector combinations TOF- 
TOF, TOF-PbSc and PbSc-PbSc, the ratio of the mea- 
sured and combinatorial like-sign invariant mass distribu- 
tions were found to be consistent with 1.0 as a function 
of the pair mass within statistical errors for all centrality 
bins. As an example, in Fig. 0] we plot the measured and 
combinatorial "+ +" and "- -" invariant mass distribu- 
tions and their ratios for the TOF-PbSc combination as 
a function of the invariant mass of K + K + and K~K~ 
pairs for minimum-bias events. As can be seen from the 
figures, these ratios are equal to the expected value of 1.0 
within the statistical fluctuations. 



The systematic uncertainty associated with this nor- 
malization procedure was estimated to vary between 0.5- 
2% for the different centralities in the different detector 
combinations. When we added all data together to derive 
the <f> spectrum, the range of the systematic uncertainty 
reduced to 0.7-1.0%. 

Finally, the cj> meson signal was obtained by subtract- 
ing the combinatorial background from the measured 
unlikc-sign invariant mass spectrum. An example of the 
K + K~ invariant mass spectrum for the TOF-PbSc com- 
bination is shown in Fig. [3] where we plotted the mea- 
sured and scaled mixed event invariant mass distributions 
for minimum-bias events. The lower panel of the figure 
shows the subtracted mass spectrum. The corrected yield 
of the 4> mesons from the experimental data is then de- 
termined by integrating the subtracted invariant mass 
spectrum within a mass window of ± 5 MeV/c 2 about 
the measured 4> mass centroid. This narrow mass window 
was used as it provided a better signal-to-background ra- 
tio compared with a wider window. Since we will show 
that there is no significant centrality dependence of the 
intrinsic width, then the extracted yields as a function of 
centrality are not being biased by the use of a constant 
integration window. The systematic effect of the mass in- 
tegration window itself on the corrected yield was studied 
by varying the stated integration limit and found to con- 
tribute 2.6-3.2%, depending on centrality, to the total 
systematic uncertainty in the integrated yield. 



E. Acceptance, Decay and Multiple Scattering 
Corrections 



The 4> meson yields were corrected for the geometrical 
acceptance of the detectors, in-flight kaon decay, multiple 
scattering effects, and nuclear interactions with materials 
in the detector using the PISA software package which is 
a GEANT-based [53| Monte Carlo detector simulation of 
the PHENIX detector. The simulation was carried out 
by generating 34 million single <p mesons in a ±0.6 rapid- 
ity interval with an exponential transverse momentum 
distribution 



dN/dpr = p T exp(~m T /(t fo + /3 2 M )) 



(1) 



with tf = 157 MeV and (3 = 0.4, i.e. an effective slope 
of T = 320 MeV in the range < p T < 10 GeV/c. The 
generated <fi mesons were then propagated through the 
simulation package. In this simulation, the BBC, DC, 
PC, TOF and PbSc detector responses were tuned to 
match the real data by including their dead areas and 
by matching their track associations and mass-squared 
resolutions. That is, the track association and mass- 
squared cut boundaries in the Monte Carlo analysis were 
parameterized to match the real data. The K + K~ pair 
acceptance efficiency as a function of transverse mass was 
calculated as 
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FIG. 4: K + K + [(a) and (a')] and K K [(b) and (b')] invariant mass spectra and ratio of real to mixed event spectra, 
respectively, for the TOF-PbSc combination. 
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The calculated acceptance efficiencies for the TOF- 
TOF, TOF-PbSc and PbSc-PbSc combinations are 
shown in Fig. as a function of tut- The points in the 
figure are located at the center of the bins. In the actual 
rriT spectra, the proper bin centroids were used. The fig- 
ure shows that the TOF detector (closed circles) has low 
acceptance for the low momentum kaon pairs due to their 



large opening angles. On the other hand, the TOF covers 
the largest mj- range for the 4> particles. As a result, the 
TOF efficiency function increases towards higher trans- 
verse mass. In contrast, the TOF-PbSc (open squares) 
and PbSc-PbSc (closed triangles) combinations offer bet- 
ter low pair momentum acceptance than the TOF. How- 
ever, the high momentum kaon identification limit in the 
PbSc leads to the efficiency function decreasing at the 
highest transverse mass values. The systematic error as- 
sociated with the acceptance correction factor originates 
from 
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FIG. 5: K + K~ invariant mass spectra for the measured 
and mixed events (top half) in the TOF-PbSc combination, 
and the subtracted mass spectrum showing the <f> meson peak 
clearly above the background (bottom half). 



i) tuning of detector alignments and mass-squared 
parameters in the Monte Carlo with reference to 
the real data (~ 3%), and 

ii) systematics in the fiducial geometries in the data 
and the Monte Carlo (~ 12%). 
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FIG. 6: Detector acceptance efficiency vs transverse mass 
of the simulated K + K~ pairs for the TOF-TOF, TOF- 
PbSc and PbSc-PbSc combinations. The statistical errors 
are smaller than the size of the data symbols. 
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F. Detector Occupancy Correction 

The high multiplicity environment in the heavy ion 
collisions produces multiple hits in a detector cell such 
as in the slats of the TOF or in the towers of the PbSc. 
These occupancy effects reduce the track reconstruction 
efficiency in central collisions compared to that in pe- 
ripheral collisions, and these occupancy dependent ef- 
fects need to accounted for in calculating the invariant 
yields. The multiplicity dependent efficiency {e occuvancy ) 
factors were calculated by embedding simulated K + K~ 
pairs into real data events. This study was done for dif- 
ferent centrality bins from to 92% in steps of 10%. 
We calculated the multiplicity dependent efficiencies for 
TOF-TOF, TOF-PbSc and PbSc-PbSc pairs separately. 
The systematic uncertainty associated with the embed- 
ding procedure was estimated for the three centrality bins 



FIG. 7: Multiplicity (occupancy) dependent efficiency cor- 
rection for detecting the <j> meson in the PHENIX detector as 
a function of the collision centrality. The most central colli- 
sions are to the right, the most peripheral collisions are to the 
left. The statistical errors are less than the size of the data 
points. 



used in the yield determinations, namely - 10%, 10 - 
40% and 40 - 92%. The systematic errors, calculated by 
estimating the occupancy efficiency corrections for differ- 
ent track confirming criteria, were found to vary from 7% 
to 10% for the three centrality bins used here, indepen- 
dent of the pair momenta. Fig. [7| shows the (-Occupancy 
factors as a function of collision centrality for the K + K~ 
pairs both identified in the TOF detector. The occu- 
pancy dependent efficiency factors were found to be in- 
dependent of the transverse momenta of the pairs. 
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IV. RESULTS AND DISCUSSION 

In this section, we present and discuss the results of 
our measurements, which consist of: 1) the (j) line shape 
analysis, 2) the transverse mass spectra analysis, 3) the 
integrated yields and ratios analysis, 4) hydrodynamical 
fits to tt^ 1 , , p, p and 4> transverse momentum spectra, 
and 5) the centrality dependence of the yields and nuclear 
modification factor Rcp as compared to that of pions and 
(anti)protons. 



Au + Au\[s^, = 200 GeV 
Same event K + K" pairs 




1.05 



1.15 



1.2 

M [GeV/c 2 ] 



A. Line Shape Analysis 

The invariant mass spectra of the <f> mesons are ob- 
tained by subtracting the combinatorial backgrounds 
from the same event K + K~ mass spectra. The details of 
the combinatorial background analysis were described in 
section III D. For the best statistical precision, we com- 
bine data from TOF and PbSc detectors to analyze the (f> 
mass centroids and widths at the five different centrality 
bins. 

Fig. |Sl shows the minimum-bias (f> — > K + K~ invariant 
mass spectrum for the PHENIX data. The subtracted 4> 
mass spectrum (lower panel), containing approximately 
5100 4> in the fit region, is fitted with a relativistic Brcit- 
Wigner (RBW) mass distribution function convolved 
with a Gaussian experimental mass resolution function. 
Using Monte Carlo studies based on the experimentally 
measured single kaon momentum resolution, the experi- 
mental 4> mass resolution is calculated to be 1.0 MeV/c 2 . 
This pair mass resolution value is found to be almost 
constant across the kinematic region of acceptance. 

The errors on the data points in Fig. [S] reflect the 
statistical errors only. The systematic errors associated 
with the mass centroid and width measurements origi- 
nate from the magnetic field uncertainties in the kaon 
momentum determination and the combinatorial back- 
ground normalization procedure. The minimum-bias line 
shape parameters (centroid and width) derived in our 
analysis are listed in Table ITT1 The fitted minimum bias 
(f> mass centroid and width are consistent at the one stan- 
dard deviation (la) level with the PDG values 2 , taking 
into account both systematic and statistical errors. 

We investigated the centrality dependence of the cf> me- 
son line shapes. For each centrality bin, we again fitted 
the cf> mass spectrum with the RBW function convolved 
with a Gaussian experimental <fi mass resolution. The re- 
sults are presented in Fig. El The left panel of the figure 
shows the centrality dependence of the fitted centroids. 



The PDG value for the <j> 
0.020 MeV/c 2 , and for the 
0.05 MeV/c 2 24] . 



mass centroid is 1019.456 ± 
d> width the value is 4.26 ± 
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FIG. 8: Minimum-bias <f) -> K + K~ invariant mass spectrum 
using the kaons identified in the PHENIX detector. The top 
panel shows the same event (circles) and combinatorial back- 
ground K + K~ mass distributions. The bottom panel shows 
the subtracted mass spectrum fitted with Relativistic Breit- 
Wigner function convolved with the Gaussian resolution func- 
tion. 



TABLE II: 4> meson mass centroid and width for the 
minimum-bias Au + Au collisions at ^/sjvjv = 200 GeV. 
The corresponding PDG values are M = 1.019456 GeV/c 2 , 
T = 4.26 MeV/c 2 



Parameter 



Value 



Centroid 
(GeV/c 2 ) 

Width 

(MeV/c 2 ) 



1.01877 ± 0.00014 (stat) ± 0.00085 (syst. 



4.24 ± 0.45 (stat) ± g;|f (syst) 



The upper and lower la systematic error limits are in- 
dicated. The dotted line shows the PDG mass centroid. 
The solid line indicates the result obtained with a one- 
parameter constant fit through the measured data points. 
These results lead to two immediate conclusions. First, 
to within less than 1 MeV/c 2 there is no observed central- 
ity dependence of the 4> meson mass centroid, and second, 
the fitted centroids at all centralities are consistent with 
the PDG value within the statistical and systematic un- 
certainties of our measurements. 

The <p mass widths, as shown in the right panel of the 
figure, are studied as a function of the centrality. The 
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error bar on each point shows the statistical error while 
the bands on the points indicate the systematic errors. 
The dotted line shows the PDG <f> mass width. The solid 
line shows the results of the constant fit assumption to 
the data points. Again, within the random and system- 
atic error limits shown, there is no convincing evidence 
of a variation of the <j) width as a function of centrality. 

The topic of medium effects on meson masses has been 
activel y in vestigated in the recent literature 0, IE El 
El SHE! Eg. The predictions are that for both cold 
and hot nuclear matter there could be a decrease of the 
<f> mass value by a few MeV/c 2 or even tens of MeV/c 2 . 
Even more dramatically the width could increase by sev- 
eral multiples above the PDG value of 4.26 MeV/c 2 . A 
sample of such predictions is given in Table IIIII How- 
ever, one of the models (2fl| considers the <fi ~> K + K~ 
channel largely insensitive to medium effects since the 
kaons are "unlikely to escape without reacting further, 
thus destroying any useful information possessed about 
the 0" . In this sense the <j> — ¥ K + K~ is inherently biased 
in that only <p decays which are unaffected by the medium 
changes, for example those produced peripherally, can be 
detected. 

It is sometimes thought that since the vacuum cr of 
the is « 45 fm/c there will be only limited sensitivity 
to medium effects in any case. However, if the resonance 
width were to actually increase by several times, as in- 
dicated in Table IIIII then the cr would then approach 
or be even substantially less than 10 fm/c, which is a 
size compatible with the expected collision volume. The 
dramatic width changes predicted in either cold or hot 
nuclear matter might be visible, at least in the dilepton 
channel if not also the K + K~ channel, as a function of 
centrality. 

The present mass centroid and width data, which are 
integrated over the available tot range, rule out any ma- 
jor changes with respect to the PDG values. Specifically, 
the one parameter fit result of 3.97 ± 0.34 MeV/c 2 ob- 
tained here excludes at the 99% confidence level a width 
value of 4.75 MeV/c 2 or greater. Possibly at the lowest 
rriT values where the (f> would remain longer in the colli- 
sion volume, or with the availability of more finely binned 
centralities, there could be visible evidence of in medium 
effects. However, the current data sample is insufficient 
to explore these possibilities. 

It seems clear from the current set of the theoretical 
models that an observed change in the <f> width would 
not be itself indicative of a QGP formation. One would 
first have to constrain the cold nuclear medium effects 
on the cj> as could be obtained in d+Au collisions, or by 
comparing peripheral Au+Au collision data results with 
central collision data. It is also important to measure 
the 4> mass in the dilepton e + e~ channel. That channel 
should be more sensitive to the <fi which are produced 
deeper or earlier in the collision volume. 



B. Spectral Shapes Analysis 

At low %, the spectral shapes carry information 
about the kinetic freeze-out conditions. Since the cen- 
trality dependence reveals the effect of the system size on 
the fireball evolution it becomes desirable to study the 
centrality dependence of the spectral shapes. Transverse 
mass spectra were obtained in three centrality bins corre- 
sponding to 0-10%, 10-40% and 40-92% of the total geo- 
metrical cross-section. We count the same event K + K~ 
pairs within a defined mass window (± 5 MeV/c 2 with re- 
spect to the 4> mass centroid) and estimate the number of 
combinatorial background pairs within that window. In 
each centrality bin, the data are divided into different mj- 
bins. The invariant mass spectrum for the same events 
and the background distributions are obtained for each 
of these tot bins. Finally, the background is subtracted 
from the same event invariant mass spectrum within the 
aforementioned ±5 MeV/c 2 (j> mass window to determine 
the number of reconstructed <j) mesons within that Tot 
bin. The reconstruction of the <f> in the Monte Carlo sim- 
ulation takes into account the effect of the <fi mass window 
size. 

The 4> mesons are reconstructed using kaons identified 
in the TOF and the PbSc detectors. Three detector com- 
binations: TOF-TOF, TOF-PbSc, and PbSc-PbSc are 
used to obtain three independent transverse mass spec- 
tra. Fig. HOI shows the minimum-bias tot spectra for the 
above three combinations. The combined result, which is 
the sum of the three combinations is also included. For 
better visibility of the data points, TOF-TOF, TOF- 
PbSc and PbSc-PbSc spectra are scaled by a factor of 
0.5, 0.1 and 0.05, respectively. The tot spectrum ob- 
tained from the PHENIX detector is fitted with the ex- 
ponential function: 

1 # N = dN / d V c -(m T -m,)/T (3) 

2-7TTOT drriTdy 2itT{T + M^) 

where dN/dy and the inverse slope T are returned as 
two fitting parameters. The lines drawn through the 
TOF-TOF, TOF-PbSc and PbSc-PbSc spectra repre- 
sent the same fit, but scaled with the same scaling fac- 
tors as the data points. Comparison of the individual 
spectra to the fit obtained from the combined spectrum 
demonstrates the consistency between the different mea- 
surements, which have different systematic uncertainties. 
The TOF-TOF, TOF-PbSc and PbSc-PbSc spectra are 
also independently fitted using Eq. (3). The resulting 
dN/dy and T are tabulated in Table llVl Both statistical 
and systematic errors are quoted. The systematic errors 
on dN / dy originate from the systematic uncertainties as- 
sociated with extraction of the yields in each tot bin (see 
Appendix C) and the uncertainties from the fitting pro- 
cedure. The latter is sensitive to the extrapolation of the 
Tot spectra to tot = rn^. A detailed account of system- 
atic errors on dN/dy and T in the full dataset from all 
sources is shown in Appendix C. We also fitted the TOF 
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TOF and TOF-PbSc data over the smaller mj- range 
of the PbSc-PbSc data and obtained consistent sets of 
dN/dy and T values from that check. 

Fig. an d Table indicate that the three differ- 
ent analyses with different systematic uncertainties give 
consistent results. This allowed us to combine the results 
and make use of the maximum available statistics in each 
rtiT bin. This combined spectrum was used to obtain the 
physics results discussed in the next sections. 

Fig. 1111 shows rriT spectra of the tfi mesons in 0-10%, 
10-40%, 40-92% and minimum bias centrality classes. 
The data points representing the invariant yield as a func- 
tion of transverse momentum are given in Appendix B. 
Each spectrum is fitted with an mr-exponential func- 
tion Eq. (3). The <fi yield per unit of rapidity (dN/dy) 
and inverse-slope (T) obtained from the fits are shown in 



Table [V] and summarized in Fig. FPU and Fig.rPol 3 . 



C. Yields and ratios 

Hadron yields and ratios carry information about the 
chemical properties of the system. The yields of strange 



3 The STAR experiment at RHIC has recently reported its analy- 
sis of the <f> K+K~ data for Au+Au at y/sjvjv = 200 GeV |6ql . 
The analyses for PHENIX and STAR have one common central- 
ity bin, 0—10%, for which the extracted dN/dy are not in agree- 
ment. The STAR value is 6.65±0.35(stat)±0.73(sj/s), compared 
with the PHENIX value 3.94 ± 0.60(stat) ± 0.62(sys). The dis- 
crepancy persists even if one eliminates the four lowest ttit data 
points from the STAR data set in order to fit over the same mj 
range for both the STAR and the PHENIX data. There is not a 
discrepancy between the quoted inverse slope parameters. 
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TABLE III: Theoretical Predictions for Changes in the <f> Resonance 



Authors 



Models 



Environment 



Predictions 



Caberara and Vacas [23] 
Pal, Ko, and Lin ^3J 



Chiral SU(3) 

AMPT 
Chiral Lagrangian 



Cold hadronic 



M4 -> 1.011 Gev 
T -> 30 MeV/c 2 



Hot hadronic M$ —* 0.95 GeV/c at twice normal nuclear density (po) 



Oset and Ramos [21] Kaon mass renormalization Cold hadronic 
Smith and Haglin [20J One boson exchange Hot hadronic 
Blaizot and Galain Nambu-Jona-Lasino Hot hadronic 



Bi and Rafelski [14] 



Bag Model 
Chiral Invariance 



Hot hadronic 



r -> 45 MeV/c 2 at 2p 
Suppression of (j) — * K + K~ relative to <f> — *• 
unchanged 
T -> 22 MeV/c 2 
M# unchanged 
-f 14-24 MeV/c 2 
M - 2Mx at T « T critical 
reduced by a factor of 6 
(/> — > K K~ disappears 
AU -> 1.029 GeV/c 2 at T « T critica i 
T -> 10 MeV/c 2 



TABLE IV: Minimum-bias dN/dy and T for different subsystem combinations. The statistical and systematic errors are shown 
after the first and second ± signs, respectively. 



Subsystem combination 


TOF-TOF 


TOF-Pbsc 


PbSc-PbSc 


Full Data Set 


dN/dy 


1.16 ±0.17 ±0.19 


1.37 ±0.15 ±0.22 


1.47 ±0.26 ±0.27 


1.34 ± 0.09 ± 0.21 


T(MeV) 


380 ± 18 ± 22 


385 ± 34 ± 28 


311 ±47 ±65 


366 ± 11 ± 18 



■ _ 

of""* 

£ 1 

0) 

S io- 1 





> 


z 

CM 

T3 


T3 

I- 

E 
















CM 



10 _ 



10" 



Au + Au @ Njs^ = 200 GeV = 

■ Full east arm - 

o TOF-TOFxO.5 

□ TOF - PbSc(east) x 0.1 

a PbSc(east) - PbSc(east) x 0.05 




0.5 



1.5 



2 



2.5 



m T - (GeV/c ) 



FIG. 10: Minimum-bias mr spectra of the measured (j) 
mesons for three different PHENIX subsystem combinations, 
with scale factors as indicated. The combined spectrum is 
fitted with an exponential function in mr, Eq. (3). The lines 
drawn through the individual spectra (TOF-TOF, TOF-PbSc 
and PbSc-PbSc) represent the same fit parameters as in the 
minimum-bias case. Statistical error bars are shown. 
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FIG. 11: m T spectra of <j> mesons for 0-10%, 10-40%, 40- 
92% and minimum-bias (0-92%) centrality classes, with scale 
factors as indicated. Each spectrum is fitted with an expo- 
nential function in mr, Eq. (3), with the fit parameters listed 
in Table V. Statistical error bars are shown. 



particles have been of particular interest as they help in 
understanding the strangeness enhancement in heavy ion 
collisions and the equilibration of strangeness. It is im- 
portant to study these phenomena as a function of system 



size. Centrality selected data can be particularly useful 
in this respect. Here we present the yield of the <j> mesons 
at mid-rapidity as a function of centrality and compare 
this yield to the yields of other hadrons and the results 
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TABLE V: dN/dy and T for different centrality bins. 



Centrality 


dN/dy 


T 


(96) 




(MeV) 


- 10% 


3.94 ± 0.60 (stat) ± 0.62 (syst) 


376 ± 24 (stat) ± 20 (syst) 


10 - 40% 


2.22 ± 0.18 (stat) ± 0.35 (syst) 


360 ± 13 (stat) ± 23 (syst) 


40 - 92% 


0.32 ± 0.03 (stat) ± 0.05 (syst) 


359 ± 15 (stat) ± 16 (syst) 


Minimum Bias 


1.34 ± 0.09 (stat) ± 0.21 (syst) 


366 ± 11 (stat) ± 18 (syst) 
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FIG. 12: Centrality dependence of 4> yield at mid-rapidity. 
The value of dN/dy increases steadily with the number of 
participants (iV par t)(left) whereas the dN/dy per participant 
pairs increases slightly from peripheral to mid-central events 
and saturates after that (right). The error bars indicate the 
statistical errors. The shaded boxes on each data point are 
the systematic errors. 



from lower energy heavy ion collisions. The dN/dy of 
(f) (shown in Fig. fTS|t is found to increase steadily with 
centrality. In the right panel, the yield is normalized 
to the number of participant pairs to take into account 
the size of the system. Within the error bars this nor- 
malized rapidity density is approximately independent 
of centrality with a possible slight increase from periph- 
eral to the mid-peripheral collisions. The trend is quite 
different from lower energy results measured at the AGS. 
In , the yield of </> was reported to be increasing faster 
than linearly with the number of participants. 

We now consider the ratio of strange to non-strange 
particles in order to understand the extent and mecha- 
nism of the strangeness enhancement in heavy ion col- 
lisions. The ratios K + /ir + and K~/ir~, are shown in 
Fig. ^] (a)-(b). Both ratios show an increase of « 60% 
from peripheral to central collisions. Most of this increase 
is for N par t < 100. Only a mild increase or saturation is 



FIG. 13: Centrality dependence of the inverse slope, T. 

observed from mid-central to the top centrality bin |5C| . 
Fig. El ( c )~(d) shows the centrality dependence of the 
0/-7T and <fi/K ratios. The limited statistics prevent us 
from extending the </> measurements into larger number 
of centrality bins spanning the more peripheral events. 
The 4>/ K ratio, in this limited number of centrality bins, 
is approximately flat as a function of centrality. The pos- 
sibility of structure in the (f>/ir ratio is difficult to infer 
from our data within the error bars. One might expect 
to see some centrality dependence in the 4>/tt ration be- 
cause there is obviously a centrality dependence to the 
K/ix ratio. However, we do not have enough centrality 
bins, nor enough signal in each bin, from which to conclu- 
sively identify a centrality dependence in 4>/i:. The flat 
behavior of the 4>/K ratio indicates that there is no pro- 
nounced difference in the production of open (kaon) and 
hidden (<fi) strangeness in heavy ion collisions at RHIC 
energies. Production in the hadronic stage via kaon coa- 
lescence K + K~ — ► <f>, seems to be excluded by that data, 
as it would result in an increase in 4>/ K~ as a function 
of N part , which is not observed. 

The collision energy dependence of the cj> yield is shown 
in Fig. 1151 where we plot the dN/dy per participant pair 
as a function of the number of participants for different 
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FIG. 14: Centrality dependence of particle ratios for (a) 
K + /tt + , (b) K~/ir~, (c) 0/0.5 (vr + + tt") (scaled by a factor 
of 5), (d) 0/0.5 (K + + K~) in Au+Au collisions at ^/s^T = 
200 GeV. 



colliding energies. The figure indicates two aspects of 
meson production at various energies. First, as we 
go from AGS to SPS to RHIC, the meson yield per 
participant increases by an order of magnitude overall. 
Secondly, at the AGS energy, we find a steady increase 
of production per participant pair from peripheral to 
central collisions. It is worth mentioning that the NA50 
experiment [2{| at CERN SPS reported an increase in 
fiducial yield (in /i + fi~ decay channel) per participant 
from peripheral to central collisions although the yield 
per participant showed saturation within error for the top 
centrality bins. The yield of mesons at RHIC, on the 
contrary, is found to be almost independent of centrality. 

To investigate further the mechanism of enhance- 
ment with increase in collision energy, we study the two 
ratios 0/7T and 0/ K~ as a function of collision energy as 
illustrated in Fig.QSI The 0/-7T ratio is found to increase 
with the collision energy from AGS to RHIC. The <j)/K~ 
ratio, on the other hand, remains almost constant within 
error bar with increasing collision energy. 
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FIG. 15: Centrality dependence of yields at different colli- 
sion energies. STAR data are from |43| . NA49 data are from 
I3H. and E917 are from El. 
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FIG. 16: 0/-7T (top half) and 4>/K~ (bottom half) ratios as a 
function of collision energy. The data symbols have the same 
meanings in both the top and the bottom halves of the figure. 



D. Hydrodynamical Model Fits to the Spectral 
Shape Data 

From the spectral data shown in the preceding sec- 
tions, we can conclude that the transverse mass distri- 
butions are well described by an exponential distribution 
and are quite similar for all the centralities. There is 
little, if any, centrality dependence of the inverse slope 
parameter in the measured centrality bins, as shown in 



Fig. El The exponential behavior is expected for particle 
production from a thermal source. 

If the system develops collective motion, particles ex- 
perience a velocity boost resulting in an additional trans- 
verse kinetic energy component. This motivates the use 
of the transverse kinetic energy, i.e. transverse mass mi- 
nus the particle rest mass, for stu dying fl ow effects. Tra- 
ditionally, the CERN experiments |57ll5S| have used sim- 
ple exponential fits to the transverse kinetic energy dis- 
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tributions and often quote just one number, the inverse 
slope T, to characterize the spectra. These fits are usu- 
ally done in the the range (tot — too) < 1 GeV/c 2 in 
order to minimize the contribution from hard processes. 

The results of such fits, obtained from previously pub- 
lished PHENIX tt ± ,K ± , P and p data Hj are shown in 
Fig. El The slope parameters show a clear mass depen- 
dence, as expected from radial flow. The mass depen- 
dence increases from peripheral to central collisions in- 
dicating stronger collectivity in the more central events. 
The <j) meson has mass similar to that of the proton. 
Hence we expect that if the 4> participates in the collec- 
tive flow then its inverse slope will be affected by this 
motion. For protons, the slope parameter changes sig- 
nificantly from peripheral to central collisions. As just 
noted, the (f> inverse slopes shows no such centrality de- 
pendence. An important difference between the results 
obtained for p, p and <fi is that the former have been fitted 
within a limited low-m-r range ( (raT^mo) < 1 GeV/c 2 ) 
as motivated above. In the case of the <f> the full measured 
range has been used for the fit due to having limited data 
at Iow-tot- As shown in Fig. ^] the three data points 
below 1 GeV/c 2 are consistent with the fit over the entire 
niT range. 
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FIG. 17: Mass and centrality dependence of inverse slope 
parameters T in mr spectra for positive (left) and negative 
(right) particles in Au+Au collisions at ^/sjvjv = 200 GeV. 
The fit ranges are 0.2 - 1.0 GeV/c 2 for pions and 0.1 - 
1.0 GeV/c 2 for kaons, protons, and anti-protons in mr — mo. 
Figure taken from reference |50||. 

The detailed study of the higher statistics (anti)proton 
spectra [5(j indicated that (anti)proton spectra cannot be 
described by a single exponential in mj — too, if the full 
measured range is considered. Although easy to visualize, 
the one-parameter inverse slope analysis proves to be too 
simplistic as a way to infer the kinetic properties of the 
system. In particular, the degree to which the </> mesons 
participate in the collective expansion will be shown (see 
Fig. 1221) to be obscured in this simple approach. 

A more sophisticated approach to this problem is to 
compare the particle spectra to a functional form which 
describes a boosted thermal source, based on relativis- 



tic hydrodynamics [59j. This is a two-parameter model, 
termed the "blast wave" model, in which the surface ra- 
dial flow velocity (/3t) and the freeze-out temperature 
{Tf ) are extracted from the invariant cross section data 
according to the equation 



dN f /p T sinhp\ /m T coshp 
- oc / f(r)rdrm T I [ — )K 1 [ 

(4) 

where Iq and K\ represent modified Bcsscl functions with 
p being the transverse boost which depends on the radial 
position according to 



p = t&ii\r 1 {(3 T )-r/R 



(5) 



Here the parameter R is the maximum radius of the ex- 
panding source at freeze-out. The function /(r) repre- 
sents the density which is taken to be uniform in this 
calculation. 

To study the parameter correlations, we make a grid 
of (Tf , /3t) pairs and then for each pair we perform 
a chi-squared minimization for each particle type. We 
use a linear velocity profile and constant particle density 
distribution. The first fit attempt is performed simulta- 
neously for the six particle species tt ± ,K ± , and the p,p 
in the range (tot — mo) < 1-0 GeV/c 2 . 

The experimental data for p and p have been corrected 
for A and A decays. However, the invariant yields of tt^ 
and K ± include feed-down from the decay of resonances 
and weak decays. To take this into account we add the 
decay of mesonic (p, r], oj, K*...) and baryonic (A, A, £...) 
resonances as follows: 

1) Generate resonances with the transverse momen- 
tum distribution determined by each combination 
of Tf and /3t- 

2) Simulate the decays using a Monte Carlo approach 
and obtain ir^ and distributions. 

3) Merge all particles, where the particle abun- 
dance is calculated with chemical parameters [6(il | 
T ch = 177 MeV, p, B = 29 MeV. 

The two-parameter Tf a vs (3t fit results obtained in 
this analysis for the most central bin are shown in Fig. 1181 
Shown in the lower panel of the figure are the x 2 contour 
levels obtained from fitting each particle spectrum sepa- 
rately. We observe that the parameters Tf Q and (3t are 
anti-correlated, the different particles have different pre- 
ferred parameter space and different sensitivity to the 
parameters. For example, the contours for the heavier 
particles are more sensitive to the flow velocity than to 
the kinetic freeze-out temperature. The minimum val- 
leys in the contours for the six particle species do overlap 
at a single common point at the 2a level. To find the 
values of the parameters at this overlap point, a simul- 
taneous fit for the six single particle spectra (n^ , , 
and p, p) was done which converges to a best fit value of 
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T fo = 108.9±|;| MeV and (3 T = 0.771±g;gg|. Using these 
parameters, we obtain the transverse momentum shapes 
shown in Fig. 1 191 where we also include the prediction for 
the 4> spectrum shape which was not part of the original 
fit. The shape of the 4> spectrum is reproduced well. 
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with the <j) data included are numerically consistent with 
the flow parameters derived previously without the <f> 
data being included. 

The two-dimensional grid search best fit values for the 
blast-wave parameterization as a function of centrality 
are tabulated in Table IVTI The radial average expansion 
velocity < (3t > is also given in this table. For the range 
of centralities studied here, the expansion velocity pa- 
rameter is seen to decrease moderately for more periph- 
eral collisions while the kinetic freeze-out temperature 
increases more significantly, approximately 40%. If one 
takes these parameters literally, then the more periph- 
eral collisions are subject to decreased radial flow while 
correspondingly the particles are decoupling kinetically 
from each other at temperatures closer to the chemical 
freeze-out temperature. This is a physically reasonable 
scenario given fewer participants in the initial expansion 
phase. 

It should be pointed out that our present <p transverse 
momentum range does not extend below 0.8 GeV/c. The 
spectral shapes at Iow-ttit, especially for the heavier par- 
ticles, are mostly sensitive to the expansion velocity. In 
the range of our <j) measurement, it is more appropriate to 
consider the asymptotic behavior of the spectral shapes, 
which for mj >> mo is given by [59|: 



T eff = 7W(l + /Jr)/(l-#r) 



(6) 



Here, T e ff is the slope parameter obtained using mr ex- 
ponential fit, as in Eq. (3). It is interesting to note that 
the measured asymptotic slopes do not seem to depend 
on centrality, although both Tf a and Pt show a clear 
centrality dependence. This is either due to a cancel- 
lation effect in Eq. since the parameters are anti- 
correlated, or indicates that the hydrodynamics descrip- 
tion is no longer valid at these large transverse momenta. 
We conclude that although the cf> data themselves can- 
not constrain the kinematic freeze-out conditions, they 
are consistent with the hydrodynamical results obtained 
from the simultaneous fit to the 7T+- , K , p and p spectra. 



FIG. 18: Contour plots for the hydrodynamical fit to the 
200 GeV single particle transverse momentum spectra for the 
K ± , and p,p in the 0-10% centrality bin. The contour 
lines are in one standard deviation steps. The upper plot is 
from a simultaneous fit with the best value shown as the dot. 
The lower plot is from independent fits for the six particle 
spectra. 

For the two other centrality bins in this study, 10-40% 
and 40-92%, we show the best fit hydrodynamical results 
in Fig. |201 and Fig. 1211 respectively. Again we see that 
the (f> transverse momentum shapes are reproduced by 
the same flow parameters which fit the identified hadron 
data at the same centrality bins. 

For a second hydrodynamical fit attempt, we include 
the <fr transverse momentum 0-10% centrality data along 
with the previously identified hadron data as part of the 
X 2 minimization search. The flow parameters derived 



TABLE VI: Blast wave model parameters 



of centrality from fitting it ,K ,p and p spectra 



as a function 
The fit 

parameters quoted here are the results from fitting the six 
identified hadrons spectra simultaneously, without including 
the d>. 



Centrality 


T f0 [MeV] 


Pt 


< Pt > 


X 2 /NDF 


0-10% 
10-40% 
40-92% 


108.9^24 
119.0tij 
150+2 


n 77 ,+U.UU3 
u - ' ' ± -0.004 
n 7/18+0-003 

0.748_ 003 

630+ 005 
u.oou_ 005 


n r 79 +U.(HJ3 
u - d ' z -0.003 

u.oou_ 002 
447+ ' 004 

U.**( -0.004 


148.0/46 
212.1/46 
150.9/46 



E. Nuclear modification factor Rcp for <j> mesons 

One of the most important result demonstrated so far 
in the study of relativistic heavy ions at RHIC is the 
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FIG. 19: Transverse momentum data and best fit hydro-dynamical results for the 0-10% centrality bin for the 200 GeV it , K , 
and p, p, along with the prediction for the <j> transverse momentum spectrum. The transverse momentum ranges for the fit are 
indicated by the solid lines, while the dashed lines indicated the extrapolated predictions for each particle species data. 



observed suppression of high p? pions in central collisions 
as compared to either pp collisions or peripheral Au + Au 
collisions 62]. A second, quite surprising observation 
is the unusually large (anti)proton-to-pion ratio at high 
pr ■ In particular PHENIX has observed [s3] that in 
central Au 4- Au collisions the p/n and p/ir ratios are 
enhanced by a factor of 3 at intermediate pr (1.5 < pr < 
4.5 GeV/c) as compared to the ratios inpp collisions and 
the ratios obtained in quark and gluon jets measured in 
e + e~ collisions [61|. It was also observed that proton 
and anti-proton production scales with N co u in this px 
region, in sharp contrast to the strong suppression of pion 
production 30] . In pp collisions high pt particles are the 
result of the fragmentation of partons. Because of the 
power law nature of the hard scattering spectrum, most 
of the particles at high pt are expected to be leading 
hadrons. The fragmentation functions, at least in the 
vacuum, are expected to be universal and independent 
of the colliding system under consideration. However, 
at intermediate p T (1.5-4.5 GeV/c) the PHENIX results 
from central Au + Au collisions are inconsistent with the 



known fragmentation functions. 

There are several conjectures which may explain the 
unexpected PHENIX result: 

1) hydrodynamic flow generated from the hadronic 
stage 65], or 

2) hydrodynamic flow generated at a partonic 
stage together with particle production from the 
recombination |32J,|33, 34] 

3) baryon junctions as a mechanism for an usually 
large build up of baryons and anti-baryons at mid- 
rapidity |63[ . 

4) in- medium modified fragmentation functions |6^| . 

The first three possibilities invoke soft processes to 
populate a region of px that is dominated by hard- 
scattering in pp collisions. The soft production at in- 
termediate pt is enhanced for protons and anti-protons, 
while pions remain dominated by hard-scattering. In 
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FIG. 20: Transverse momentum data and best fit hydrodynamical results for the 10-40% centrality bin for the 200 GeV 
7v^~ , K ± , and p,p, along with the prediction for the cj> transverse momentum spectrum. 



conjecture (4), the production for both pions and p,p is 
by hard-scattering, but the fragmentation functions are 
modified in the medium, thus resulting in an enhanced 
p/t: and p/it ratios. 

The second of these conjectures is particularly impor- 
tant since this hypothesis presupposes a partonic state 
with collective behavior. A critical factor which may 
differentiate among these scenarios is whether the large 
abundance of protons is due to its mass, or to the number 
of constituent valence quarks |22, |33|, OS |33 . Hydrody- 
namic flow generated at the hadronic stage imparts a sin- 
gle velocity to the moving matter, hence similar mass par- 
ticles should exhibit the same momentum increase from 
this effect. In contrast, hypotheses 2 and 3 are dependent 
on whether the particles are baryons or mesons. 

The (f> meson which has a mass similar to that of a 
proton yet, like the pion, has two valence quarks, should 
distinguish between (1) and (2) or (3). We examine 
the scaling properties of the intermediate px yields of 
the <j> an d compare those to the yields of (p + p)/2. 
Fig. shows the transverse momentum spectra mea- 



sured in three different centralities, each scaled down by 
the number of binary nucleon-nucleon collisions. The 
(anti)proton spectra show two pronounced features. Be- 
low pt < 1-5 GeV/c, the spectral shapes are strongly 
influenced by the radial flow and thus the more central 
data have harder slope. Above pt — 1-5 GeV/c, the ef- 
fect of radial flow is negligible. The spectra converge to 
the same line. Moreover, they scale with N co u for all cen- 
trality classes, as expected for hard-scattering unaffected 
by the nuclear medium. The <$> spectra have a quite dif- 
ferent behavior. There is no visible curvature at lower pt 
since this curvature is not expected to be prominent in 
the measured region. At higher px the <f> spectra run par- 
allel to the (anti)proton spectra, but do not obey N co u 
scaling. To examine this feature on a linear scale, we 
plot the ratio between the central and peripheral data, 
i.e. the ratio Rcp (Fig. I23|) . The systematic error aris- 
ing from the determination of N co u is represented by the 
dotted bar and is about 19%. This systematic error is 
common for all three particle species shown in the figure. 
The solid bar around Rcp = 1 represents the N co u error 
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FIG. 21: Transverse momentum data and best fit hydrodynamical results for the 40-92% centrality bin for the 200 GeV 
7v^~ , K , and p,p, along with the prediction for the cj> transverse momentum spectrum. 



for the protons (the same 19%). Since these errors are 
correlated, if the curves in the figure were to change due 
to a change in the value of N co n, they will move together 
within the extent of the error bars shown. We would like 
to emphasize the comparison between the Rc p values for 
the protons and the <fi. In this comparison, the system- 
atic errors in determining N co ii cancel. The important 
systematic errors to consider are those that can move 
the 4> points with respect to the proton points. When 
determining a ratio of spectra measured at different cen- 
tralizes, most systematic errors cancel. After removing 
the N co a error, the sources of error that remain for the 
(j> come from the multiplicity dependent corrections and 
the effect of the mass window as described above. For the 
protons and pions |29l l30| , the error that remains that is 
independent from the error determined for the </>, comes 
from the multiplicity dependent corrections. The rela- 
tive error between the <j> and the proton measurements 
is evaluated at 7% and is represented by the extended 
solid bar just below Rcp = 1. Clearly the <j>'s behavior 
is more like that of the pions and not like that of the 



protons. Thus we conclude that the <fi meson exhibits 
a suppression effect at intermediate pr similar to that 
of the pions. Although we can not conclude whether <j> 
production at this intermediate pr is dominated by soft 
or hard processes, this observation provides support for 
models which depend on the number of valence quarks in 
the particle as opposed to models which depend upon just 
the mass of the particle in order to explain the anomalous 
proton yields. 



V. CONCLUSIONS 

A systematic measurement has been made of 4> 
production at mid-rapidity in Au+Au collisions with 
yfsWN = 200 GeV at RHIC. The general fea- 
tures of the data include the yield which rises from 
0.318±0.028(stat)±0.051(syst) in peripheral collisions to 
3.94±0.60(stat)±0.62(syst) in central collisions. There 
is seen to be little centrality dependence to the inverse 
slope which is about 360 MeV. The centroid mass and res- 
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FIG. 22: pt spectra of proton and <j) mesons at different 
centralities scaled down by their respective number of N co u 



onance width are extracted with high enough precision 
to rule out any large (MeV/c 2 scale) deviations from the 
accepted PDG values. At pr below 1.5 GeV/c, a blast 
wave analysis of the most central pion, kaon and proton 
spectra with a freeze-out temperature Tf D of 109 MeV, 
and a transverse velocity (3t of 0.77 describes the most 
central 4> data as well. A similar conclusion, with mono- 
tonically changing results for Tf Q and (3t, holds for the 
less central events. At higher transverse momenta, all 
particles lie above the blast wave fits, which suggests that 
the dominant particle production mechanism is no longer 
soft physics but is giving way to the expected jet frag- 
mentation. A study of the ratio Rcp of the (f> provides 
a critical new piece of information in understanding the 
anomalous proton-to-pion ratio seen in central heavy ion 
collisions at RHIC, since the is a meson with a mass 
similar to that of a proton. The Rcp value for the <j) 
above a transverse momentum of 1.5 GeV/c is about 0.6, 
similar to that of the pions but inconsistent with the pro- 
ton value of 1. This indicates that the </> meson is being 
suppressed in this pt range for the more central colli- 
sions. The lower px blast wave fits imply hydrodynamic 
behavior at the hadronic stage, while the inconsistency 
with simple hydrodynamics at higher pt, is something 
one would expect in the jet fragmentation region. This 
transition is an important factor to consider as one be- 
gins to understand the mechanism of particle production 
in central collisions at RHIC. 
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FIG. 23: Neon scaled central to peripheral ratio Rcp for (p + p)/2, ty° , and <j>. The proton and pion results are published |3(1 |. 
The vertical dotted bar on the right represents the error on 

N coii° % / N coU 92% ■ The shaded solid bar around RCP = 1 represents 
12% systematic error which can move the proton and/or <j> points with respect to one another. The dotted horizontal line at 
Rcp = 0.62 is a straight line fit to the <j) data. 
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The number of - - pairs is: 
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We therefore get: 
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(N(n) - (n 2 )) (A6) 



The last expression is trivial: the number of like-sign 
pairs is equal to the total number pairs minus the number 
of of unlike-sign pairs. 

The ratio R of unlike to like pairs is therefore: 



R = 



n+. 



N(n) 



n ++ + n. 



N(N — l)/2 — (N(n) - (n 2 )) 



(A7) 

This is a general result, free of any assumption and should 
always be fulfilled, in pure combinatorial background as 
well as in a mixture of signal + combinatorial back- 
ground. 



APPENDIX A: COMBINATORIAL 
BACKGROUND 



2. Combinatorial background 



Assume there are N tracks per event of which n are 
positive and N-n are negative, the probability of the par- 
tition being P(n) such that J2n=a P( n ) = 1- 

All expressions below refer to the average number of 
pairs per event 



Consider a pure combinatorial background sample. By 
its essence there are no correlations between tracks in 
such a sample i.e. within an event the probability p to 
have a positive (or a negative) track is constant and in- 
dependent of the number of tracks. Therefore the prob- 
ability P(n) of having n positive tracks out of the total 
of Af tracks is given by a binomial distribution: 



General relation between the number of like 
and unlike pairs 



P(n) 



n\{N-n)\ 



p n (l-p) 



N-n 



(A8) 



The total number of pairs is: 

N . _ N(N~l) 

~ * 'nn.'/.rs 



(Al) 



(n) = pN 



(A9) 



(n z ) = o z + (n) = Np{\ -p)+ p 2 N 2 (AlO) 
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Replacing these values in expression (6) gives: 

ip-p 2 ) 



R 



(1/2 -p + p 2 ) 



(All) 



If there is charge symmetry, i.e. p = 0.5 one gets R = 1 
and consequently the combinatorial background is given 
by: 



n.?^ = n++ + Ti- 



(A12) 



This is an exact relation. It holds with quite good accu- 
racy even if there is some charge asymmetry. For example 
for an asymmetry of 10% (i.e. p = 0.525)R=0.995. 



3. The formula N?? = 2y/N++N- 



The combinatorial background is rigorously given by 
this formula provided that the number N of tracks per 
event has a Poisson distribution: 



V(N) 



(N) 



N„-N 



M 



(A13) 



Again we assume that the N tracks are divided into n 
positive and (N-n) negative tracks, the partition is given 
by the binomial distribution (A7) and all expressions be- 
low refer to average number of pairs per event. 
Number of + + pairs: 



oo N , -i \ 

AT ++ =£p(A0 5>(n)^P (A14) 

N=2 n=0 



Using relations (A7) and (All) and some algebra leads 
to: 



N ++ = ~i>W 
Similarly the number of - - pairs is given by: 

oo N 



tv__ = 5>(ao]>>W 

N=2 n=0 



(N — n)(N — n - 1) 



-(i-p)W 



(A15) 
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The number of combinatorial background + - pairs is: 



JV 



N=2 n=0 



= p(l-p)(N 2 ) 
Inspecting (A15), (A16) and (A17) shows that: 



N 



CD 



= 2y/N ++ N. 



(A18) 



(A19) 



(A20) 



APPENDIX B: DATA TABLES OF CENTRALITY 
SELECTED A SPECTRA 



The invariant yields, 



of the <b mesons 



2irmT ' dniTdy 1 

in different centrality bins are shown in Table fVTTl 



APPENDIX C: SYSTEMATIC ERRORS ON dN/dy 
AND T 



The sources of the systematic errors on yield (dN/dy 
and T) measurements are as follows: 

a) Systematic error on the combinatorial background 
normalization, S norm : This originates from the sys- 
tematics of the event mixing. Since both same event 
and mixed event like sign distributions represent pure 
combinatorials, we estimated unlike sign combinato- 
rial background by normalizing the mixed event unlike 

sign distributions to 2,1 Nf^ event ■ event and 



2^JN^ xed event ■ N^i xed event and the difference in the 

extracted <f> signal from the real data between these two 
normalizations are attributed as the systematic uncer- 
tainty. 

We use the same normalization factor for all mj- bins. 
So, the above systematics are applicable to the measured 
dN/dy only, not on the inverse slope, T. 

b) Systematics of <j) mass window, S mass : We count 
the number of reconstructed <p mesons by integrating the 
(j) meson invariant mass spectra within ± 5 MeV mass 
window with respect to the measured centroids in both 
data and Monte Carlo. The systematic associated with 
this mass window is estimated by measuring the extent of 
the changes in dN/dy and T after constructing <fi meson 
m,T spectra within five different mass windows ± 3, ± 5, 
± 8, ± 10 and ± 15 MeV, with respect to the measured 
4> centroids. 

c) Uncertainties in extrapolation of <f> meson mr spec- 
tra to hit = W0 ; S ex trap- This is studied by 

(i) fitting the mj- distributions using exponential and 
Boltzmann functions, two different fitting functions, ex- 
ponential and Boltzmann, and 

(ii) fitting the transverse mass spectra within different 
hit ranges. These are applied to both dN/dy and T. 

d) Acceptance correction systematics, 5mc'- The sys- 
tematics associated with acceptance correction factors 
derived from Monte Carlo analysis are investigated by 
considering two sources: 

(i) Tuning of detector alignments in Monte Carlo with 
reference to the real data (~ 3%), and 

(ii) Systematics in the fiducial geometries in data and 
Monte Carlo (~ 12%). 

The systematic error from this source is independent 
of the momenta of reconstructed <fi mesons. So, this is 
attributed to dN/dy only. 
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TABLE VII: m,T spectra of cj> mesons in different centrality bins. The systematic errors on invariant yields are from combina- 
torial background normalization, <j> counting mass window, acceptance correction efficiencies from Monte Carlo and occupancy 
dependent corrections. 



Centrality 


rriT 


m,T bin size 


2 

1 d N 
2-KrriT dmj'dy 

[(GeV 2 /c 4 )- 2 ] 


Stat, error 


Syst. error 


(%) 


(GeV/c 2 ) 


(GeV/c 2 ) 


[(GeV 2 /c 4 )- 2 ] 


[(GeV 2 /c 4 )- 2 ] 


- 10 


1.365 


0.4 


0.51114 


0.15120 


0.0778747 




1.691 


0.2 


0.19114 


0.04656 


0.029121 




1.891 


0.2 


0.12065 


0.02702 


0.018382 




2.091 


0.2 


0.07313 


0.01758 


0.011142 




2.291 


0.2 


0.03335 


0.01372 


0.005081 




2.565 


0.4 


0.01683 


0.00541 


0.002564 




2.891 


0.2 


0.01353 


0.00406 


0.002061 




3.300 


1.0 


0.00268 


0.00074 


0.000408 


10 - 40 


1.364 


0.4 


0.37232 


0.04145 


0.05667 




1.691 


0.2 


0.09224 


0.01255 


0.01404 




1.891 


0.2 


0.04872 


0.00719 


0.00742 




2.091 


0.2 


0.03867 


0.00483 


0.00589 




2.291 


0.2 


0.02320 


0.00389 


0.00353 




2.564 


0.4 


0.01246 


0.00160 


0.00189 




2.891 


0.2 


0.00415 


0.00110 


0.00063 




3.294 


1.0 


001 20 


00021 


00018 


40 - 92 


1.364 


0.4 


0.04432 


0.00622 


0.00645 




1.691 


0.2 


0.01295 


0.00185 


0.00189 




1.891 


0.2 


0.01009 


0.00110 


0.00147 




2.091 


0.2 


0.00531 


0.00074 


0.00077 




2.291 


0.2 


0.00287 


0.00060 


0.00042 




2.563 


0.4 


0.00126 


0.00024 


0.00018 




2.891 


0.2 


0.00062 


0.00019 


0.00009 




3.293 


1.0 


0.00019 


0.00004 


0.00003 



TABLE VIII: Systematic error in dN/dy. 



Centrality 

(%) 


$ norm 

(%) 


&rnass 

(%) 


&extrap 

(%) 


Smc 
(%) 


$occu 

(%) 


riot 
®sys 

(%) 


Minimum Bias 


0.7 


2.6 


4.0 


12.4 


7.8 


15 


- 10 


0.8 


4.2 


5.2 


12.4 


10 


17 


10 - 40 


1.1 


2.3 


5.9 


12.4 


8.5 


16 


40 - 92 


0.7 


3.1 


6.0 


12.4 


7 


16 



TABLE IX: Systematic error in T. 



Centrality 


&mass 


Sfit 


fitot 
u sys 


(%) 


(%) 


(%) 


(%) 


Minimum Bias 


0.6 


4.9 


5 


- 10 


1.1 


5.2 


5 


10 - 40 


1.1 


6.2 


6 


40 - 92 


1.1 


4.2 


4 



e) Systematic error in the occupancy dependent effi- 
ciency corrections, S occu : The systematic error associated 
with this efficiency is estimated by calculating the occu- 
pancy dependent correction with different track confir- 
mation criteria and is independent of the pair momenta. 
This is a systematic effect on dN/dy only. 

The above systematic errors are quoted in Tables lyTiTl 
and llXI for dN/dy and T, respectively. 
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